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CLASSICAL RELAXATION OF BRg MOLECULES IN AN ARGON HEAT BATH 

by Ronald Razner 
Lewis Research Center 

SUMMARY 

Monte Carlo computer calculations of three-dimensional classical trajectories were 
carried out for collisions of bromine (Brg) molecules with argon (Ar) atoms in an argon 
heat bath at 1800 K. Limited calculations were also performed at heat bath temperatures 
between 300 and 10 000 K. The results were used to study energy transfer mechanisms 
and to obtain probability distributions for internal energy and angular momentum changes 
per collision. These indicate that large changes in the internal energy of Brg are pre- 
dominantly accompanied by large transfers in the angular momentum of the molecule. 
Starting with Brg molecules at zero angular momentum in a heat bath above 1000 K, a 
significant fraction of these molecules will acquire in a single collision an angular 
momentum large enough so that vibration- rotation coupling in the Brg molecule cannot 
be ignored. This implies that molecules with a significant vibration- rotation coupling 
will be common collision partners with argon atoms in successive collisions. This fur- 
ther implies that simple mechanisms of pure translation-vibration and pure translation- 
rotation energy transfer which are used in many theories of rate processes may not be 
adequate for the B^-Ar system. The relevance of these results to shock tube and other 
kinetic studies is discussed. 


INTRODUCTION 

Recent shock tube (refs. 1 to 3) and flash photolysis (refs. 4 to 6) experiments in- 
volving diatomic molecules in a heat bath of inert gas atoms have led to discussions re- 
garding the validity of several theoretical models for chemical rate processes (refs. 1, 
5, 7, and 8). The object of this work was to elucidate possible energy transfer mech- 
anisms that may have been neglected in simplified theories of molecular relaxation. 

For this purpose classical three-dimensional collision trajectories were calculated for 
the bromine-argon (Brg-Ar) system, and probability distributions for internal energy 
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changes and angular momentum transfers were evaluated. These probability distribu- 
tions and also individual trajectories were studied to determine important features of 
such collisions. In particular, the role of angular momentum and vibration- rotation 
coupling in the relaxation of the bromine molecule was investigated in some detail. 

The calculations were done specifically for the case of Brg molecules in the ground 
electronic state embedded in a heat bath of Ar atoms, where the population of Brg mole- 
cules is small enough for Brg-Brg encounters to be neglected. There are three main 
reasons for choosing this system for investigation: (1) current interest generated by 
availability of extensive shock tube and flash photolysis data, (2) recent spectroscopic 
work on Brg (ref. 9) which has led to the determination of an improved ground state in- 
teraction potential (ref. 10), and (3) the fact that the ground state of Brg consists of a 
large number of closely spaced vibrational levels (ref. 11). The last feature is de- 
sirable from the standpoint of classical mechanics, where the energy levels for a mole- 
cule are taken to be continuous. Thus, one might expect that such system would provide 
an optimum chance for classical mechanics being relevant to the physical situation. 

Throughout this report statistical tabulations of pertinent results are presented in 
considerable detail, since these can be of use in attempts to formulate analytical models 
for energy transfer processes. Unless otherwise stated, all data are presented in 
scaled SI units defined as follows: 


Unit 

Symbol 

Value in SI units 

Time 

t 

10" 14 sec 

Length 

o 

A 

10 meters 

Mass 

m 

10" 26 kg 

Velocity 

V 

10 meters/sec 

Energy 

u 

10" 18 joules 

17 
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(1 joule = 10 ergs 


= [6. 241(97)±(12)]10 18 eV 
- [5. 034(80)±(18)]10 22 cm -1 
= [1. 4399(8)±(4)]10 20 kcal/mole) 


Temperature 

T K 

Angular momentum 

Jm | 2 (with or m 2 A^v 2 

squared 

without 


subscripts) 
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COLLISION MODEL AND CALCULATION METHODS 


Figure 1 shows the three-body coordinate system, where B-C represents the 
diatomic molecule and A the atom. The total Hamiltonian governing the three-body 
motion is 

H = H CM + H ABC 

where 



with reduced masses given by p B( -, = m B m c /m BC and p^ BC = m A m Bc/ m ABC* The 
Hcm * s the kinetic energy of the A -B-C center of mass and does not contribute to the 
energy exchange in a collision. The vector Rg^, locates particle B relative to par- 
ticle C, while the vector X locates the atom relative to the center of mass of B-C. 

The center of mass of Br 0 was taken to be the geometrical center of the molecule. 
Actually, a natural sample of bromine consists of Br (50. 54 percent) and Br 
(49.46 percent), so that the center of mass of Brg corresponds to the geometrical center 
for only about half the molecules. It is doubtful that the error is significant here be- 
cause of the large mass of Brg compared to mass differences between isotopic species. 
However, the reduced masses p B( ~, and p^ B C were calculated from the average 

masses of natural samples of the elements (79.909 and 39.948 atomic mass units for 

12 

Br and Ar, respectively, on the C scale). 

Interactions are taken to act along lines joining the particles and are given by 
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v 2 (i 5 acI)= 4e [(yi5 A ci) 12 - (vi5 AC i) 6 ] 


(5) 


Equation (3) is a potential of the form suggested in reference 12 which was fitted to an 
RKR potential for the Br-Br interaction by the method outlined in appendix B. Equa- 
tions (4) and (5) representing the atom-molecule interaction are Lennard- Jones (6-12) 
potentials with constants e and a estimated from diffusion data (ref. 13) by mixing 
rules (private communication with Roger A. Svehla, NASA Lewis Research Center). 

Hamilton’s equations for the relative A-B-C motion are given by 


R BC " 9H / SP BC ~ P Bc/^BC 

-P BC = 9H / 3R bC = 1 BC aV l(l^Bc|)/ 3 | R Bcl " \ ^B 0 V 2 (l^AB l) I R A b| 




(6) 


x = 3H/ap = p/p A)BC 

-p = 8H/SX =T ab 3V 2 (|r ab |)/3|r ab | +T ac 3v 2 (|r ac |)/s|r ac | ^ 


where Pgc an ^ P are momenta conjugate to R B( ~, and X, respectively, and 
‘jk'V'V are unit vectors. With e = 0 the equations would describe the motion of 
a free atom and a free homonuclear vibrating rotator with any vibration- rotation inter- 
action that may be present. 

For given initial conditions the equations of motion (in Cartesian coordinates) were 
integrated on an IBM 7094 computer usually employing a variable step-size routine de- 
veloped by Dugan, Rice and Magee (ref. 14). This routine uses the principle of con- 
servation of energy to choose an optimum step-size at each time-step interval 

in the numerical integration. (Conservation of the square of angular momentum about 
the A-B-C center of mass could also be used. ) One decides beforehand on a maximum 
variation in H^ B £ (in this case 0. 1 percent) that can be tolerated for a complete tra- 
jectory. The routine then breaks this up into step tolerances, taking large step sizes 
in noncritical prats of the integration, and smaller step sizes at critical parts. If after 
any individual step the incremental error is greater than it should be at that point, the 
calculation for that step is rejected and repeated with the step size halved. On the other 
hand, if the incremental error is smaller than the tolerance for a given step, succeeding 
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steps are chosen larger. In actual practice high energy conditions were more easily 
handled by the computer than low energy conditions (that is, a smaller portion of the 
allowable error was used up). Overall, this routine has proven to be much more effi- 
cient than using a fixed step size. For certain conditions to be subsequently considered, 
trajectories could not be calculated with sufficient accuracy by using the variable step- 
size routine. In such cases a small fixed time-step interval of 0.25 t was used. The 
correctness of the two integration methods was verified by checking time reversal in- 
variance. 

The allowed tolerance on H^g^ determined noise levels for the integration; that 
is, it established lower limits for changes in a molecule's internal energy and angular 
momentum. Changes larger than these lower limits could be considered significant, 
while changes smaller than these lower limits could be regarded as zero. By inspection 
of individual trajectories it was found that noise levels of 0. 0001 u for energy and 0. 001 

2 o 2 2 

m A v for angular momentum squared could be conservatively set for all trajectories 
studied. 

An interaction radius for the atom-molecule interaction (see fig. 1) was determined 
empirically by calculating a series of trajectories with the atom placed initially at vari- 
ous distances from the B-C center of mass, and the molecule given a series of specified 
initial orientations. The initial translational motion of the atom was always along the 
vector -X for these trajectories. The initial value finally chosen for |x| was 10 A. 

For larger initial values of |x| no changes in trajectories with the same initial B-C 
orientation could be detected above the noise level. In subsequent calculations integra- 
tion was carried out from |x | = 10 A (with |x | decreasing in time) to |x J = 10 A (with jx | 
increasing in time). The latter defined final conditions for a trajectory. For each tra- 
jectory the following quantities were recorded: 


(1) H^BC’ final, (f), and initial/final (i/f) 

(2) Square of total angular momentum about A-B-C center of mass, 

where M = (X x P) + (Rg^ x (f) and (i/f) 



= M . M 


(3) Impact parameter b defined as the perpendicular distance of the atom trajectory 


from the B-C center of mass that would result if e = 0 


(4) |X| ; (f) and (i/f) ^ ^ ^ ^ 

(5) Scattering angle for atom (0 to v), cos~*[x(i)/ X(f)/|x(i)| • |x(f)|] 

(6) Total internal energy of B-C, Hg C = M BC |Rg C | 2 + v i(I^bcO’ and 

(7) Vibrational energy of B-C, H y = ( 2 M BC )~ * P Bc) 2 + V lO R BcD’ 

and (i/f) 

(8) Rotational energy of B-C, H r = I m bC^ ^^BC^^Bcl ^ where M BC = R BC XP BC’ 

(f) and (i/f) __ __ 

(9) Square of angular momentum of B-C, ^ = Mg^, • Mg^,; (f) and (i/f) 
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(10) Kinetic energy of relative motion, M A BC ' X I ’ (f) and (i/f) 

(11) Angle between initial and final B-C angular momentum vectors when 

2 /J ' "-" irr ' r M„o(f)/|M nn (i)| * |M RC (f)|] 


| Mg(-, | ^(i) =£0, cos jMg£i(i) * .M.gQvi.'/ i x, ' l BC VJ ' / 1 ' i' l “BC'‘ 

Initial conditions were either assigned values or else were randomly selected from 
an appropriate distribution function by means of the formulas derived in appendix C. 

The molecule was assigned an initial energy Hg^, and also a magnitude for its angular 
momentum vector. The orientation of the angular momentum vector and the orientation 
of the molecule were then randomly distributed over a sphere. Internal energy and 
angular momentum were given assigned initial vaiues in preference to vibrational energy 
H y and rotational energy H r because the former are constants of the motion for the 
free molecule, while the latter are not. 

With the atom initially at a distance of 10 A from the center of mass of the molecule 
(see fig. 1), the magnitude of the relative momentum |p| was chosen from a Maxwell - 
Boltzmann flux distribution corresponding to a specified heat bath temperature T, while 
the orientation of this momentum was randomly distributed over a half sphere. The 
impact parameter b was calculated from the initial conditions and compared to a pre- 
assigned maximum impact parameter b . When b ^ b the trajectory was cal- 

niax max 

culated and the collision termed a real collision (RCOL). When b > b the trajec- 

max 

tory was not calculated and the collision was called a null collision (NCOL). These 
latter cases, although not calculated, were nevertheless counted as events for the de- 
termination of the distributions for energy and angular momentum changes in a collision. 

The maximum impact parameter was calculated as b_„„ = 4 liL l + 2 1/6 a. 


where R 


l BC 


b — = -!| r B( j + 2 1 / 6 


max 


max 


is the outer classical turning point in the B-C vibration and the 


max 

second term locates the minimum in the Br-Br interaction. 


This choice of b was 
max 


made because trajectories with b > b were found, by trial, not to contribute sig- 

max 

nificantly to the overall energy transfer. 

A statistical sample of initial atom energies H^(i) must be distributed as a 
Maxwell-Boltzmann flux distribution (ref. 15) 


dw = 



for RCOL alone, for NCOL alone, and for RCOL and NCOL taken together. Table I 
shows how well this was accomplished for a typical set of data consisting of 5229 total 
collisions, of which 1756 were RCOL (corresponding to the row labeled Hg^-,(i)/D e = 0. 5, 
| M bc | 2 (i) = 0 in table III). 
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II 



In scaled units the physical constants relevant to this report are as follows (with 
i = (-1) 1 / 2 ): 

k = 0. 138054xl0' 4 uK' 1 

K 2 = 0. 111197xl0' 3 u 2 t 2 
D x = 0.31801 u 

D 2 = 0. 57547x10' 3 - iO. 88651xl0~ 3 u 
Dg = 0. 57547x10' 3 + iO. 88651xl0~ 3 u 
o? 1 = 0. 19756X10 1 A' 1 

n 2 = 0. 15994X10 1 + iO. 64990X10 1 A' 1 

a 3 = 0.15994X10 1 - i0.64990xl0 1 A ' 1 

(Z = D e = 0. 31916 u 

a = 0.3867X10 1 A 
e = 0. 3411x10' 2 u 
|R BC (e)| = 0. 22816 lxlO 1 A 

M b c = 0.663420X10 1 m 

Pa bc = 0. 530667x10* m 
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VIBRATION- ROTATION COUPLING IN THE FREE Br 2 MOLECULE 

As observed in the preceding section, the rotational energy H r and vibrational 

energy H y are not constants of the motion for the free diatomic molecule, whereas 

the internal energy H or , = H + H is a constant of the motion. This is due to the fact 

BL v r 

that H r is proportional to 1/1 Rgc ' > which varies with the vibrational motion. Asa 
result, whenever the molecule possesses a finite angular momentum, a periodic coupling 
between H y and H r exists in phase with the vibrational motion. Thus, as a prelimi- 
nary step to considering atom-molecule collisions, the amplitude of vibration- rotation 
coupling was calculated for a range of internal energies and angular momenta and com- 
pared to quantum mechanical vibrational level spacings (J = 0) obtained from RKR cal- 
culations (see appendix B). It is likely that the phenomenon of vibration- rotation 
coupling, shown schematically in figure 2, may become significant for an atom-diatomic 
collision theory whenever the coupling amplitude in the free molecule is comparable to 
energy spacings in successive RKR levels. 

For a measure of the H y , H r coupling we can use the peak-to-peak amplitude of 
vibrational energy H y (p - p) defined by 

H V (P - P) = H y> max - H v> min 


where H 


is the instantaneous vibrational energy corresponding to the outer 


classical turning point in the B-C vibration, and H 


corresponds to the inner clas- 


sical turning point. This amplitude was determined at constant H B( ~, as a function of 
| | ^ for the free Br 2 molecule from dummy trajectories obtained by setting e = 0 

in the equations of motion and conserving Hg^, to a tolerance of 0.05 percent. The 
results are shown in figure 3 for Hg^, = 0.1 D g , 0. 5 D g , and 0. 9 D 0 . (The curve 
labeled H Br in fig. 3 is considered separately at the end of this section. ) One 

Jt5 \ ? IHaA 2 

point together with its coordinates is shown on each curve. From 


M. 


M. 


BC 


BC i =0 up to 
2 to within the 


this point H y (p - p)/H B( -. increases linearly with increasing 
limits set by the noise level in the integration. 

A list of RKR (J = 0) vibrational levels (see appendix B) for the ground state of the 

7 Q Ol 

Br Br molecule is given in table II. Using this table together with figure 3, one can 
determine the extent of H y , H r coupling and can estimate when a theory of energy 
transfer conceived as a process of transitions between pure vibrational levels is likely 
to fail. For example, a coupling amplitude of H y (p - p)/H B( ~, = 0.01 corresponds ap- 
proximately to 1/20 of a RKR vibrational level spacing at HgQ = 0. 1 D g , 1/3 of a level 
spacing at 0. 5 D , and 1 level spacing at 0.9 D . 
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Average values for H v and H r were also determined for the free molecule (at 
constant H B( -,) as a function of |m B( -,| 2 by integrating the vibration and rotation energy 
functions according to the relation 


<H 



(7) 


At constant H B( ~,, 


<*V " h bc (I^bcI / 2m bc) 4^1 ^bc I 2 !) 


(8) 


where <j R BC | 2 )> is the mean square of the internuclear separation. To a good approxi- 
mation the numerical data can be expressed by a linear relation: 


<n v > " h bc x I m bcI 


0) 


where x = 0. 0143 


at Hgc = 0. 1 D g , x = 0. 0134 


at H BC " °* 5 D e’ 


and x = 0. 0084 at 


Hbc = 0.9 D g (with energies and angular momentum expressed in scaled units). Cor- 
responding constant values of ^iRgC^ are 5.27(1), 5.62(5), and 8.9(7) K-^, respec- 


(Uncertain figures are enclosed in parenthesis. ) The linear relation holds for 

i ^ i 2 

Ht,^ = 0.1 D„ and 0. 5 D as far as the calculations were taken (up to | Mbc I = 10 

e ■- ■ 2 


tively. 

[ BC 

or H y = 0), whereas for H B( -, = 0.9 D g the linear relation is valid up to iMgcl* = 1.0. 
For larger angular momenta at Hg^-, = 0.9 D g the numerical value of x increases 
slowly to 0. 01 and then varies irregularly between 0. 0099 < x < 0. 0130 with increasing 
IM-ooI 2 . This increase of x with Im^^I 2 reflects the fact that at constant an 


BC 
increase of 


BC 1 


M t 


and hence a decrease in H 


BC 
The 


iiit/i caot; I J-vx-r-j I produces an increase in H, auu a u m xx . xn& 

oO /, — r . 2\ v 

smaller H v , however, implies a smaller <J R B q I /• The irregularities in x at large 

| Mgc | ^ are perhaps an indication that some errors are being accumulated in the inte- 
gration of the equations of motion at high rotational energies. However, in all cases 
the average rotational and vibrational energies satisfy Hg^, = (H v ) + <H r >. 

Also shown in figure 3 is a H v , H r coupling curve for the molecule near disso- 
ciation, which for present purposes is defined to be H — 0.9999 D P . (It is not 

v , max .S. . 

practical to define dissociation as H ^ D because the condition ]R T -. r ,| — °° can- 

V 2 niaA c DL 

not be handled by the computer. ) With this definition the maximum internal energy 
H^p, that the molecule can have without dissociating is determined from dummy 

dL> y m cLx 

trajectories to be 
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( 10 ) 


H BC,ma>e < 0 - 9999 + y|«Bcl 2 

where energies and angular momentum squared are expressed in scaled units and 
y = 0. 0100±0. 0001(m 2 A 2 v 2 ) - The H . HL coupling curve for is almost 

V x Dv y iiloA 

linear, in contrast to the other curves in figure 3, because the former was calculated 

for constant H instead of for constant H Dr ,. 

v, max jdu 

SIGNIFICANT FEATURES OF INDIVIDUAL Br 2 -Ar COLLISIONS 

A typical atom- molecule trajectory for close encounters is shown in figure 4, where 
energies and distances are plotted as functions of time. (This particular example was 
taken from data with initial conditions Hg^.(i) = 0.1 D e , |Mg£ | 2 (i) = 0, and T = 1800 K. ) 
Most of the energy transfer occurred within a small time interval, comparable in mag- 
nitude to the time required for one cycle in the Br 2 vibration. In terms of the atom- 
molecule interaction this time interval corresponds to Br-Ar distances where V 2 be- 
comes repulsive. (One cannot conclude from this that the attractive terms in V 2 may 
be ignored. The statistical influence of these terms is considered later in this report. ) 
Also evident in the figure is the existence of a significant H v ,H r coupling, which was 
discussed in the preceding section. In this section the influence of atom velocity on 
impulsive collisions and the sensitivity of trajectories to changes in poetntial constants 
are considered. 


Dependence of Impulsive Collisions on Initial Atom Energy 

Many simple atom-diatomic collision theories make use of an impulsive collision 
approximation, where atom- to- molecule energy transfer is considered to take place in 
the duration of one molecular vibration period or less. For the present collision model 
we have investigated whether a collision that is impulsive at high heat bath temperatures 
would remain impulsive as the heat bath temperature is lowered. For this purpose a 
series of trajectories with identical impact parameters and Br 0 initial conditions was 
calculated by scaling the components of initial atom velocity, Xj, X 2 , Xg (see appen- 
dix C) for the trajectory in figure 4 to correspond to heat bath temperatures from 
300 to 10 000 K where 


X.(T) = (T/1800) 1 / 2 X.(1800) 


(ID 
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Results are shown for T = 300, 1000, 1800, and 10 000 K in figures 5(a) to (d), respec- 
tively, where the variables ,H BC’ and H r are P lotted the same way as in fig- 

ure 4. Since the initial frequency of the B-C vibration (as indicated by |R B q| ) is the 
same in all these figures, one sees that if a collision is impulsive at T > 1000 K, it 
will likely remain impulsive as the atom velocity is decreased to correspond to a heat 
bath temperature as low as 300 K. However, there are limitations to this statement 
that are considered in the concluding discussion of this report. 


Sensitivity of Trajectories to Potential Constants 

An attempt was made to estimate the sensitivity of individual trajectories to changes 
in potential constants by recalculating the trajectory shown in figure 4 under the following 
conditions: (1) use of a simple Morse potential instead of equation (3) where Dj = D g , 

D 2 = Dg = 0, ttj = 0. 19689X10 1 A -1 ; (2) e = 0. 2411x10" 2 u; (3) e = 0. 4411X10" 2 u; 

(4) a = 0.4367x10* A; and (5) cr = 0.3367x10* A. In each case, all constants not listed 
had the original values. These changes in the atom-molecule interaction V 2 are prob- 
ably more drastic than any changes that could be expected from improved experimental 
data. For example, even for other inert gases we estimate the constants: 
e = 0. 995x10" ^ u, a = 0. 3438x10* A for He; e = 0. 1878xl0" 2 u, a = 0. 3544x10* A for Ne; 
e = 0. 4072xl0~ 2 u, a = 0. 3984x10* A for Kr; e = 0. 4597x10” 2 u, a = 0. 4180x10* A for 
Xe; and almost all of these lie within the previous variations in e and a. Final values 
and the ratio (i/f) are listed for all pertinent quantities - first for the original trajectory 
in figure 4 then for cases (1) to (5): 

Impact parameter: same for all trajectories; b = 0.2338711 A 

Scattering angle (radians): 2.6092621, (1) 2.6104993, (2) 2.6244507, (3) 2.5905407, 

(4) 2.6970121, (5) 2.4579556 

**ABC^ : same for a ** trajectories; H AB( -,(f) = 0.09076595 u (original) 

H ABC (i/f): t- 0000893 * (1) 1.0001728, (2) 1.0001322, (3) 1.0001293, (4) 1.0001100, 

(5) 1.0002308 

H BC (i/ f ): 0.5682207, (1) 0.5626052, (2) 0.5719623, (3) 0.5632753, (4) 0.6167731, 

(5) 0. 5173918 

H^(i): same for all trajectories; H A (f) = 0.03459789 u (original) 

H A (i/f): 1.7012077, (1) 1.7501041, (2) 1.6830058, (3) 1.7265491, (4) 1.5097518, 

(5) 2.0231369 

|m bc | 2 (f): 1.6896222 m 2 A 2 v, (1) 1.6551220, (2) 1.5350854, (3) 1.8447367, 

(4) 1.7241471, (5) 1.7011569 

The use of the Morse potential resulted in displacing both the original maximum 
(0. 0357 u) and minimum (0. 0286 u) in the final H y upward by 0. 0019 u. 
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In this report, the concern is mainly with the total energy transferred to the mole- 

i -► | 9 

cule AHg£ and also the angular momentum transfer A | M B( -, | . These are considered 

in the next section. To significantly alter any statistical results to be discussed, AHt-,/-. 

■ — 1 2 

and A| Mg£ | revisions due to changes in potential constants of at least a half order of 
magnitude for individual trajectories would have to add up in a noncompensating way for 
a major portion of the real collisions (RCOL) involved in a random sampling. I believe 
it is safe to say that any such revision is extremely unlikely, and hence that the constants 
e and a are adequately determined for present purposes. 


Estimate of Three-Body Interactions 

The interactions (3) to (5) chosen for the present collision model assume only a 
pairwise additivity of forces acting along lines joining particles A, B, and C. In dense 
systems such as compressed gases, liquids, and solids, pairwise interactions may no 
longer be adequate, and a consideration of three-body and many-body forces might then 
be necessary. Hence, the question may arise as to what influence a three-body potential 
term, Vg(|R^ B |, |r B q|, I r ^q|)j would have on individual trajectories. 

To estimate the possible influence of three-body forces a potential of the form of a 
van der Waals interaction between three atoms was used (ref. 16): 

V„ = (1+3 cos A cos B cos C) (12) 

(PabI • I 5 bcI • PacI ) 3 

where cos A, cos B, and cos C are cosines of the angles BAC, ABC, and ACB, re- 
spectively, and the constant v was assumed to have the value v = 7. 32 uA which holds 
for Ar-Ar-Ar interactions (ref. 17). For the trajectory in figure 4, Vg came to within 
0. 1 percent of the total energy (0. 001 H^ B( -,) for only a few time units during the closest 
approach of the atom to the molecule. At the same time, however, the interaction Vg 
was more than 150 times that of Vg. Consequently, it appears that if the value used 
for v is at least of the correct order of magnitude, then no significant error will result 
by completely neglecting three- body forces. 
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CALCULATIONS FOR Br 2 MOLECULES IN AN ARGON HEAT BATH 

Probability Distributions for AHg^, and 

Probability distributions for the molecule undergoing a change AHg^, a| Mg^, | ^ 
per collision with an argon atom randomly selected from the heat bath were determined 
for given initial values of internal energy Hg^(i) and angular momentum squared 
|MbcI^U)' A family of such distributions is illustrated schematically in figure 6. 

Each distribution curve depends on all the initial conditions 

f < iH BC ) = f ( AH BC; H BC (1) - I“bc| 2(1) ’ t ) 


t ( a l 5l Bc| 2 )= t ( A l i5l Bc| 2 ; H BC<»’ K C | 2 <‘>> T ) 

so that each set of initial conditions ^Hg£,(i), |Mg^,|^(i), t) determines a specific distri- 
bution for AHg£ and a specific distribution for A|Mg£,| From a data set consisting 
of 5229 total collisions, of which 1756 were RCOL, calculated for initial conditions 
(0. 5 D , 0, 1800 K), it was found that data sets consisting of 1000 RCOL were sufficient 
to determine any given distribution. For an initial angular momentum |M^^| (i) = 0, 
the distribution ^ would consist of only a "gains" branch. A statistical 

equilibrium for or A | | ^ would result if "gains" and "losses" branches for 

the corresponding distributions were symmetric, whereas a relaxation would occur if 
either a "gains" or a "losses" branch were predominant. Several types of equilibrium 
conditions which are relevant to subsequent discussions can be defined as follows, where 
n stands for noise level: 

(1) Overall gains-losses balance for AH 


/ -n /~+ 0 ° 

„ f < AH BC> d(AH BC> = _/ +n f < AH BC> d < AH BC> 


( 13 ) 


— 2 

(2) Overall gains-losses balance for A Mg^, 



(3) Simultaneous gains-losses balance for AHg^, and a|m B q| 2 when conditions 

(1) and (2) are both satisfied 

(4) Statistical equilibrium for AHg^, (gains and losses branches symmetric) 



) = 



f(AH BC )d(AH BC ) 


for all a s n 


(5) Statistical equilibrium for a| Mg^ 


(15) 



/ a + ” f ( i Kc| 2 K A l“Bc| 2 ) 


for all a ^ n (16) 


(6) Simultaneous statistical equilibrium for AHg^ and A|Mg£,| 2 when conditions 
(4) and (5) are both satisfied 

Tables III and TV list probability distributions for several combinations of initial 
conditions in terms of the fraction of total trajectories resulting in changes in a given 
AH B c or a I m bC^ range* From these tables several results are apparent. For 
most initial conditions, a major portion of the total trajectories is within the noise level 
of the calculation, and relatively few trajectories are responsible for the bulk of the 
energy and angular momentum transfer. Referring to table III as an example, one finds 
that collisions of the latter type are distributed so that, with |m b ^| 2 (i) = 0 and Hg^i) 
from 0 to 0.9 D g , after a single collision 28 to 38 percent of the molecules have 
| Mg£ | 2 5: 0. 001, 23 to 29 percent | Mg^ | 2 ^ 0. 01, 17 to 21 percent | Mg^ | 2 2: 0. 1, 
and 4 to 9 percent | M B< ~, j 2 ^ 1. 0. In terms of the quantum formula | Mg^, | 2 = R^J(J + 1) 
these correspond (to the nearest integer) to transitions from J = 0 to J 2: 3, 9, 30, 
and 95, respectively. 

It should be emphasized that all fractional populations in tables III and IV are given 
on a total collision basis. These would change if the radius of the sphere of interaction 
were changed, because increasing the radius of the sphere of interaction would increase 
the proportion of NCOL. On the other hand, the relative distribution of RCOL in given 
AHg£ and A | Mg C | 2 ranges does not depend on the radius of interaction provided this 
is chosen large enough so that the atom and molecule can be considered initially free. 

By using the ratios of real collisions to total collisions in the tables, the fractional 
populations as given can be converted to a RCOL basis. Thus, for example, if 5.8 per- 
cent of the molecules with Hg^(i) = 0. 5 D g and | Mg^ | 2 (i) =0 at T = 1800 K have 
| Mg C | 2 £ 1.0 after a single collision, this represents 17. 3 percent of the real colli- 
sions . 
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Calculations where the internal energy of the molecule H B( ~, is close to the disso- 
ciation energy D g present a special problem in that the molecule can dissociate during 
a trajectory. Thus, final conditions for such trajectories must be defined as |x| = 10 A 
with |x| increasing in time, or H y ^ 0.9999 D g , whichever occurs first (see eq. (10)). 
Trajectories of this type were calculated for T = 1800 K, and the resulting probability 
distributions are included in table III. For each initial angular momentum (i), 

initial energies were set at H B ^-,(i) = 0.999 Hg^ max (see eq. (10)), and trajectories 
were calculated using a constant time-step size to avoid false dissociations by accumu- 
lation of integration errors. For these calculations b m was set at 7.966 A, which 

corresponds to 2 1 / 6 ct + 4 - (Ror'l for H = 0.9999 D . This value of b mnv , 

z | £> '- / 1 max v , max e max 

was found by trial to include all those trajectories capable of supplying the small incre- 
ment of energy needed for dissociation. From table HI it is seen that the probability 
of a molecule dissociating increases significantly as |m b ^-,| ^(i) increases from zero. 

The nature of angular momentum changes in the dissociation process is shown in fig- 
ure 7, which compares the frequency of events where dissociation proceeds with an in- 
crease in |M Rr ,| with the frequency of dissociations accompained by a decrease in 
1 -*■ | 2 

|M B cl • The table in the figure shows the noise level population of dissociation events. 
For small initial angular momenta dissociation events with a gain in angular momentum 
outnumber dissociation events with a loss in angular momentum. As | M B( -, | z (i) increases 
the loss-to-gain ratio for dissociation events reaches a maximum, and then decreases. 


Sensitivity of AH B( ~, and a|m B q|^ Distributions to Heat Bath Temperature 

When the results in table IV are plotted, the behaviors of AHgg and a|m b( -,| % 
probability distributions as functions of heat bath temperature can be obtained for initial 
conditions ^H B( -,(i), | M B( ~, | ^(i)J = (0, 0), (0. 5 D 0 , 0), and (0. 5 D g , 1. 0). Such plots are 
shown in figures 8(a) and (b). It appears that distributions over a wide range of heat 
bath temperatures can be obtained with reasonable confidence by interpolation once cal- 
culations have been done for a few temperatures. It is also seen that for AH., ,, and 
,o 

A|M bc | ranges which include small to intermediate changes, there is a significant 
difference in the temperature dependence of probability distributions, depending on 
whether or not the molecule has an initial angular momentum. The most pronounced 
differences are in the heat bath temperature range of 300 to 2000 K. In this region 
AHg£ and A|M b ^-,| gains for initial conditions (0,0) and (0. 5 D g , 0) are decreasing 
functions of temperature, while those for (0. 5 D g , 1. 0) are increasing functions of tem- 
perature. The opposite dependence is found for AHg^, losses for (0. 5 D 0 0) and 
(0. 5 D g , 1. 0). Thus, in the temperature range of greatest interest (corresponding to 
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most of the available experimental data) it would appear that a consideration of angular 
momentum transfers as well as internal energy changes is important. 


Relation Between 


AH BC 


and A | M 



Distributions 


For given initial conditions the probability distributions f(AHg^) and f [A | | j 

are not independent of each other; that is, given a specific energy change AHg^, all 
A | | ^ changes in the distribution f^A|M B Q | ^ are not allowed and vice versa. The 

correspondence between internal energy and angular momentum changes was determined 
for each data set in tables III and IV. Results are presented as gain-loss matrices in 
tables V, VI, and VII, where entries represent fractions of RCOL (including changes 
within the noise level). 

Table V correlates overall gains and losses for T = 1800 K, where it is seen that 
diagonal changes ^AH B( -,, A | | ^ = (g, g) and (l, l) are the most probable. Much less 

probable are off-diagonal changes, the least probable being those in the (g, l) quadrant 
(with the possible exception of conditions near dissociation which are complicated by 
dissociation events). Diagonal events appear to undergo a smooth transition with in- 
creasing | MgQ | ^(i) from a dominance in the (g, g) quadrant to a dominance in the (l, l ) 
quadrant. Also, as | M B( ~, | ^(i) increases from zero, events in the (l, g) quadrant appear 
to undergo a steady decrease (conditions near dissociation excluded). However, the 
least probable off-diagonal quadrant (g, l) behaves differently from the other three in 

I -+■ 1 O 

that its population first increases with increasing |M B( -,| (i), reaches a maximum, and 
then decreases. Finally, with the possible exception of conditions near dissociation, an 
equilibrium or balance in diagonal events does not coincide with a balance in off-diagonal 
events. The latter, if it occurs (not certain with present calculations), occurs at a 
higher |m B( -,| ^(i) than the balance for diagonal events. 

Table VI consists of a breakdown of each quadrant in the matrices of table V into 
subquadrants which separate large changes from small changes. Table VII is a similar 
compilation for the data sets in table IV. The dividing line for energy was set at 
AHbc = 0. 01 u, which is about three times the zero point energy. (Vibrational transi- 
tions from RKR level 0 to RKR levels 1 and 2 would correspond to AH y = 0. 0054 u and 
0. 0127 u, respectively. ) Internal energy changes less than this were arbitrarily called 
small, and those at or above 0.01 u were called large. Similarly, the dividing line for 
A I | was set at 0.1 m A v (or J = 30). Partitions labeled 1 correspond to small 
energy and small angular momentum changes; those labeled 2 correspond to small 
energy and large angular momentum changes; those labeled 3 correspond to large energy 
and small angular momentum changes; and those labeled 4 correspond to large energy 
and large angular momentum changes. Partitions 1 include all those changes which 
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correspond to | Aq| <1 for H , where q is the RKR vibrational quantum number, 
whereas partitions 3 would include all pure (or almost pure) vibrational energy changes 
where I Aql >2. For all conditions it is seen that the probability of the latter is almost 
negligible. Of all RCOL, the most probable changes are large and small energy gains 
(losses) accompained by significant angular momentum gains (losses), as represented 
by the sum of partitions 2 ^ and 4^ for angular momentum gains and the sum of parti- 
tions 2^ and 4^ for angular momentum losses. 


Equilibrium as a Function of Heat Bath Temperature 


Sets of initial conditions, (HjjqU), | | ^(i), x) , 


were determined which result in 
2 


overall gains-losses balances for AH B( -, and/or A|M B( -,|“ by graphically locating 
solutions to equations (13) and (14) which define these balances. Results are presented 
in table VIII for Hg^(i) = 0. 5 D g through a heat bath temperature range of 300 to 
10 000 K, and for H B( -,(i) = 0. 1 D g and 0. 9 D g at a heat bath temperature of 1800 K. 

Also given in table VIII is a comparison of the vibrating- rotator balance points with mean 
angular momenta for rigid rotators at various intemuclear spacings according to equi- 
librium statistical mechanics. Since I R B c I ^ ^ or a rotator is constant, the mean 

internal energy is simply kT per degree of freedom, or <H r ) = kT = | ^}/ 

2/igc |r bc | 2. From this table it appears that if one were given initial conditions 

HptrJi) and T for the vibrating- rotator and wished to estimate a corresponding value 

r 1 * 1 2 | — 1 2 

of M Dr , (i) that would result in a A M nr , balance, the required estimate for 
oC £><~ . 2\ 

T < 1800 K would correspond moderately well to <\ |M B £,| >for a rigid rotator with an 

internuclear spacing equal to the average internuclear spacing for the vibrating- rotator 

at the given energy. (See discussion accompanying eq. (9).) For example, a a|m b( -,|^ 

balance occurs at | M 00 | ^(i) = 0. 98 for the vibrating rotator with H xir ,(i) = 0.5D„ and 

dc 1 , o\ e 

T = 1000 K. The corresponding ( \ Mg^, | for the rigid rotator is 1.03, a difference of 

only 5 percent. However, as the temperature increases beyond 1800 K, the vibrating 
rotator balance occurs at angular momenta corresponding to a rigid rotator with de- 
creasing internuclear spacing. At T = 10 000 K this spacing is approximately equal to 
the inner classical turning point for the vibrating rotator. 

Whereas the equilibrium angular momentum for a rigid rotator is simply a constant 
times the heat bath temperature, one finds that for the vibrating rotator a plot of values 
of |M B £,|^(i), which yield a|m b ^,|^ balances as functions of heat bath temperature, 
has a curvature that becomes more pronounced as the heat bath temperature increases. 

If at constant H B( -,(i) a similar plot is constructed for values of |Mg£,| (i) correspond- 
ing to AHgc balances, the AHg^-, and a|m B( -,|^ curves intersect at a heat bath tem- 
perature corresponding to a simultaneous gains-losses balance for AHg^, and 
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a|M;bc| 2 * Interestingly enough, this temperature corresponds approximately to a 
statistical equilibrium temperature for internal energy obtained by setting 

H BC«) = I kT. 

Investigations similar to the previous one could perhaps be carried out for internal 
energies near dissociation, but here the situation is complicated by the fact that mole- 
cules are dissociating with a probability of dissociation that depends on | Mg^, | (i). 
Calculations performed so far are not extensive enough to consider such conditions fur- 
ther in this report. 

The information presented in this and preceding sections permits an estimate of 
the magnitude of H y , H r coupling likely to exist in the Brg molecule as it collides with 
an Ar atom. Even if a collection of molecules with zero initial angular momentum were 
embedded in a heat bath at T > 1000 K, after a few collisions a considerable fraction of 
these molecules would possess angular momenta of the order |Mg£ | > 1. If the angu- 

lar momenta corresponding to gains-losses balances in table VIII are regarded as typical 
for Brg molecules under the specified conditions, then from figure 3 it is seen that a 
H y , H r coupling amplitude of the order of 10 percent or more of Hg^, would be com- 
mon. At internal energies Hg^ > 0. 5 D g this corresponds to a coupling amplitude 
well in excess of the RKR vibrational level spacing. 


Scattering Angle and Reorientation of Mg^, 

The atom scattering angle and also the angle between the initial and final B-C angu- 

i ^ i 2 

lar momentum vectors (for |Mg£,| / 0) were tabulated and studied. Unfortunately I 

did not decide to record the latter quantity until this work was well under way, and hence 
results cannot be given for all data sets in tables III and IV. However, sufficient re- 
sults are included for a meaningful study of trends. 

For all data sets, the scattering angle populations decreased monotonically with in- 
creasing scattering angle. As might be expected, the decrease was most rapid at the 
highest heat bath temperatures, showing that atom trajectories for high energy collisions 
are effectively straight lines. 

For planar motion (two-dimensional theories) the Mg^ reorientation angle would 
be restricted to the values 0 or it, whereas in the present three-dimensional collision 
model this angle can have any value in the interval [0,7 r]. Tabulations of the Mg^ re- 
orientation angle together with AHg^, and a| Mg^, | 2 indicate that a large angle almost 
always implies large energy and angular momentum changes. Events consisting of a 
large reorientation of the Mg^ vector with little or no net change in Hg^ and |Mg^-,| 
were extremely rare. 
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The way in which the Mg^, reorientation angle populations were distributed (for 
RCOL) is shown by bar graphs in figure 9 for all data sets where this quantity was re- 
corded. Each bar graph indicates the angle intervals representing approximately 50, 

70, and 90 percent of the data. (The exact fractions as tabulated are given above each 
bar graph, since tabulations were only performed for steps of 0. 1 radian. ) Some trends 
are discernible for data taken at fixed Hg^/i) and | Mg^, | ^(i) through a range of tem- 
peratures and at fixed Hg^(i) and temperature through a range of | Mg^, | ^(i). Data for 
(0. 5 D g , 1.0, T) show a slight increase in the range of Mg^, reorientation angles with 
increasing temperature up to 3000 K, but a restriction to angles of 0. 1 radian or less 
at 10 000 K. This restriction to small angles at 10 000 K holds for all angular momenta 
that were considered, except for a sharp increase for the data set (0. 5 D , 8. 24, 

10000 K), which is a gains-losses balance point for AlMg^l . This may indicate the 
possibility of a resonance at gains-losses balance conditions, where the probability of 
large MgQ reorientation angles increases. Such possibility for AHg^ and/or 
A | Mg£ | 2 balance conditions is suggested by comparing the data sets (0. 5 D e , 0. 25, 

300 K), (0.5D o , 0.82, 1000 K), (0.5D o , 1.48, 1800 K), and (0.9 D 0 , 2.0, 1800 K) 
with others in figure 9. Finally, the data sets near dissociation at 1800 K show a de- 
crease in the range of Mg^, reorientation angles as the angular momentum increases. 
(The three bar graphs in fig. 9 labeled "attractive term in V £ removed" are dis- 
cussed below. ) 

Plots of scattering angle as a function of the Mg^ reorientation angle showed a 
random scatter of data points, and hence no apparent statistical correlation between 
these two quantities. 


Effect of Removing the Attractive Part of Vg 

Some theories neglect the attractive term in a potential of the form V 2 . Hence, I 
have evaluated the effect of the attractive terms by removing -(cr/|R^g|)6 anc * 

- (a / 1 R^q |)^ from the atom- molecule interaction and recalculating three probability 
distributions - (0.5 D„, 0.25, 300 K), (0.5D fl , 1.48, 1800 K) and (0. 5 D , 8.24, 

10 000 K). Results are presented in table IX in the form of probability distributions and 
in the form of matrix tabulations as previously discussed. Since e/k = 247 K, one 
might expect the greatest influence to be for heat bath temperatures close to this value. 
A comparison of table IX with corresponding entries in tables IV and VII confirms this 
expectation. In particular, the calculations at 300 K show almost a perfect gains- 
losses balance for A|Mg£ | in contrast to a predominance in gains for calculations 
with the attractive Vg terms included. At 1800 K what was previously a gains-losses 
balance for AHg^, is now converted to a dominance in losses. Results for 10 000 K 
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are unchanged. Finally, figure 9 includes tabulations for the Mg^-, reorientation angle, 
which again shows a large influence of the attractive interaction at 300 K. 


RESULTS AND DISCUSSION 


Diatomic molecule-third body collision theories of energy transfer are often con- 
ceived in terms of translation- rotation processes (rotational relaxation) and translation- 
vibration processes (vibrational relaxation) where the two pure transfer modes are as- 
sumed to be independent to a good approximation. (A recent review of vibrational and 
rotational relaxation processes is given in ref. 18.) The former process is character- 
istic of rigid rotator theories, and the latter process of collinear collision models. 

Most of the theories of dissociation of diatomics involve in some way the idea of ex- 
citing an unperturbed oscillator, such as the quantum mechanical process of vibrational 
ladder climbing, with an activation energy for dissociation equated to the dissociation 
energy. However, attempts to correlate theory with experiment (ref. 1) have indicated 
that an activation energy is probably less than the dissociation energy. This being the 
case, energy transfer mechanisms other than a pure translation-vibration process must 
be sought. 

From an inspection of individual trajectories (e.g. , see fig. 5) we have concluded 
that one may, to a good approximation, assume an impulsive energy transfer between 
the atom and the diatomic molecule. Such an assumption is common to many theories 
and has been used, for example, in the quantum mechanical calculations of Bak and 
Fisher (ref. 19). However, this does not necessarily mean that such impulsive energy 
transfer can be construed as a pure vibrational or pure rotational excitation. Rather, 
it appears that a sizeable portion of the trajectories responsible for the bulk of the 
energy transfer (RCOL) will involve molecules with a significant vibration- rotation 
coupling, commonly exceeding one or more RKR vibrational level spacings. As a con- 
sequence, the molecules would be so perturbed that one could not justify considering a 
dissociation process in terms of a pure translation- vibration energy exchange. For 
example, it has been pointed out by Bak and Fisher (ref. 19) that if activation and dis- 
sociation energies are not equal, one could not assume a random walk on the energy 
levels of an unperturbed oscillator. At best the process would have to be considered as 


a H a to H bc 


transfer with vibrational and rotational modes in the B-C motion cross- 


feeding each other. Dissociation would presumably have to be defined as H — D , 

° v, max e 

but at a given average vibrational energy (H ) sufficient angular momentum could be 


added to the molecule to make it dissociate at the outer classical turning point in the 


vibrational motion. In addition, a comparison of the way probability distributions for 
|m bc I ^(i) = 0 and | Mg^ | ^(i) » 0 vary with heat bath temperature indicates signifi- 
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cant differences. Hence, it would not be surprising for theories that ignore the influence 
of angular momentum to yield significantly different dissociation rates, compared to 
theories that include the effect of angular momentum. Any discrepancies between vari- 
ous theories should increase as the average number of collisions per dissociation event 
increases, because each collision is effectively a 
governed by distributions f(AHg ( -,) and f^A|Mg^ 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 9, 1970, 
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SYMBOLS 

scaled unit of length, 10” meters 
speed of light 

dissociation energy for B^ 

potential constants for Br - Br infraction, £d. = D 

J e 

Hamiltonian (total energy) 
kinetic energy of relative Brg-Ar motion 
energy of relative motion of particles A, B, and C 
internal energy of Brg 

kinetic energy of translation of Brg-Ar center of mass 
rotational energy 
vibrational energy 

2 

Planck’s constant squared divided by 4n 
unit vector 

rotational quantum number 
Boltzmann's constant 

angular momentum about A-B-C center of mass 
angular momentum about Brg center of mass 

scaled unit of mass, 10 kg, subscripts indicate masses of particles 

relative Brg-Ar momentum 

momentum of particle B relative to particle C 

vibrational quantum number 

location of particle B relative to particle C 

equilibrium internuclear separation for Brg 

temperature, 

- 14 

scaled unit of time, 10 sec 

- 18 

scaled unit of energy, 10 joules 
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V(x) 

v^IRbcD 

v 2 (|r A bI),v 2 <I 5 acI> 

v 3 (|5 A bU5 bc I.I5 ac I) 

V 

X 

Z 


“m 

e 

0 

6 

M 

cr 

v 


<P 

X^, w \T , 0)1 
o’ o o’ o J o’ el 

co x , co y , w z„ | 
e e’ e J e’ e e J 


interaction potential 

interaction potential for Br-Br 

atom- molecule interaction potentials 

three- body potential 

scaled unit of velocity, 10^ m/sec 

location of Ar relative to Brg center of mass 

symbol for a random number in the interval [0, 1] (with 
or without subscripts) 

potential constants in Br-Br interaction, Vj 
potential constant in Morse potential 
Lennard-Jones potential constant in Vg 
longitude in cylindrical coordinates 
colatitude in spherical polar coordinates 
reduced mass (with or without subscripts) 
Lennard-Jones potential constant in Vg 
potential constant for Vg 
longitude in spherical polar coordinates 

molecular constants 
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APPENDIX B 


CURVE FITTING PROCEDURE FOR V^R^I) 

The potential V^Rg^l ) (see eq. (3)) is the generalization of the Morse function 



suggested in reference 12 that can be made to coincide with "true" RKR potential calcula- 
tions for diatomic molecules to a degree not possible with a simple Morse function. 

This generalized potential is of the form of a sum of three Morse functions 



where x= |Rg£,| - |Rgc( e )| anc * D g =X)Dj is the dissociation energy referred to 

V(0). In equation (B2) pairs of D^ and pairs of Oj may occur as complex conjugates, 

with the restriction Re(cVj) > 0 to ensure that the potential remains finite as x — °°. 

In attempting to construct a potential of the form of equation (B2) we have found that 

literally following the procedure given in reference 12 will not unequivocally lead to 

79 8 1 

physically acceptable solutions. For ’ Brg both older (ref. 9) and recently revised 

molecular constants (ref. 10) led to complex a. with negative real parts. The situation 

79 79 81 81 J 

was not changed for ’ B^ and ’ Brg where molecular constants were estimated 
from the isotope effect (ref. 11). However, given RKR calculations, one can replace 
the first part of the calculation in reference 12 with a direct method of curve fitting. 

This yields a similar set of final equations for the D^ and where certain numerical 
criteria established in reference 12 (to be explained) are not required. Using the spec- 
troscopic constants of LeRoy and Burns (ref. 10), we have determined V(x) correspond- 

79 81 

ing to the ground state of ’ Brg. The RKR potential was first recalculated with a 
program available at this laboratory (ref. 20). The value of D (16068 cm - *) was ob- 

— I ^ 

tained from the experimental value for D of 15906 cm (ref. 9) plus the zero point 
- 1 ° 

energy of 162 cm . Our calculations agree closely with those in reference 10, except 
in the following particulars. There is a very slight discrepancy in that a calculation of 
spectroscopic energy as a function of vibrational quantum number (for J = 0) using the 
constants given in reference 10 possesses a maximum at 16053. 7 cm'\ which is in- 
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consistent with the dissociation energy of 16068 cm’^. Our RKR calculations give 
78 vibrational levels up to the energy 16053. 7 cm'* (see table II) instead of the 80 levels 
reported by LeRoy and Burns. Also we find |Rg^-,(e)| = 2. 28161 A compared to LeRoy 
and Burns’ value of 2. 2814 A. / -a x \2 

Let V(x q ) be the potential energy at a given RKR level q and Fj(Xq) = (l-e •* . 

Then ' ' 


D l F l (x q ) + D 2 F 2 (x q ) + D 3 F 3 (x q } = V(x q } 


(B3) 


where F^(x q ) can be expanded as 

F i ( V%^[ <2n - 2)/n: ](-“fq) n 


Substituting (B4) into (B3) and retaining terms up to x 
D l[(' a l X q) 2 + (” a l X q) 3 + ( 7 / 12 )(-«i x q ) 4 + 


6 

q 


f 


gives 

+ (62/720)(-a lXq ) 6 ] 


(B4) 


+ D 2 [(-“ 2 x ,) 2 + — ] + D 3 [(“ 3 x q ) 2 + — ] “ V( V <B5> 

Defining G^ = + BgOfg, with k = 0, 2, 3, 4, 5, 6, Gq = D e , and the other 

G^ obtained from RKR data, equation (B5) becomes 


( x q) 2 G 2 - ( x q) 3 °3 + < 7 / 12 >( x q) 4 G 4 ‘ W(\T G 6 + (62/720)(x q ) 6 Gf . = V(x q ) 

(B6) 


where « has been replaced by = for the purpose of curve fitting. Five values of x q (for 
classical turning points) and corresponding values of V(x q ) are now chosen so as to be 
representative of the RKR calculations over the range of energies for which an accurate 
curve fit is desired. This gives five equations which are then treated as simultaneous 
equations and solved for numerical values of G 2 to Gg. These values of G^ plus Gq 
are then used to determine the D. and a. from solutions to the dimensionless equa- 
tions: 

(k= 0,2, 3, 4, 5, 6) (B7) 

J=1 
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where bp = 1, = (3^/0^^, ' z ] ~ D j/ D e ’ an< * = a j^ a M‘ These equations are 

solved for and if ^ by the method suggested in reference 12. The value of is 
calculated from 


°M = K + w o x o ) ( 7TC/i/RD e) 1/2 (B8) 

where c is the speed of light and j u the reduced mass. 

In the calculations of reference 12, equation (B8) is the first step in an iterative 
procedure for a(^a^) from u> Q + u) Q x 0 + pco Q y 0 = a^hD e / 7 TCjuy*/ 2 where p is an iter- 
ation parameter and neglecting pa> Q y o gives a = a^. For the case of Br 2 the iteration 
converges with the initial step for a considerable numerical range of p. This may be a 
source of difficulty. The quantities oj q , W 0 X 0 > are calculated from the molecular 

constants co e , u Q x e , u e y g , o> e z e by 

“o = w e ‘ w e x e + ( 3 / 4 ) w e y e + (1 / 2)w e z e 

w o x o = w e x e “ ( 3 / 2 ) co e y e " (3 / 2 ) w e z e 

“o y o = ^eJ^e + 2w e z e 

For 8 *Br 9 we have (in cm - *), o> = 324. 24, w x = 1. 172, w y = 0.00342, and 
a> e z e = -0.000101 (ref. 10). 

For present purposes we have chosen two RKR points for |Rg(-.| < |Rg£,(e)| an< ^ 
three points for |Rg^| > l^gc( e )| (Q = 18_, 2_, 0 + , 4 + , 18 + ) for which values of b^ were 
determined as 

b Q =l b 4 = 0.6743075 

b 2 = 0. 9960699 bg = 3. 441186 

b 3 = 0. 7493079 bg = 8. 840423 

In general there will be two sets of solutions for and <f., but in the present instance 
only one set was physically acceptable. In the procedure of reference 12 bg = b 2 = 1 
always, b 2 = b 4 , and bg and bg are sought as close to unity as possible, such crite- 
rion being in principle necessary for convergence in these authors’ original expansion. 

In our case we obtain values for b 2 , bg, and b 4 close to those obtained by the method 
of reference 12 (although b 2 / b 4 ); but, values for bg and bg considerably exceed 
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unity. The D. and corresponding to these are given in the text (see p. 7). The 
potential so determined was found to be a significant improvement over the Morse po- 
tential for low vibrational energies, but no improvement above q = 18. This curve fit 
was a compromise between a desire to achieve as close a fit as possible for low lying 
RKR vibrational levels and a desire to represent the required curvature on both sides of 
|R BC ( e )| f° r high vibrational levels. By fitting the potential to more points for any de- 
sired portion of the RKR calculations one can, of course, achieve a better fit for that 
portion. For example, choosing q= 10_, 0 + , 3 + , 5 + , 10 + resulted in a fit for 
I r bcI > jRBc( e )l to within 1 percent of the RKR calculations up to q = 60, but a poor 
fit for |R BC | < |R BC (e)|. 
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APPENDIX C 


SELECTION OF INITIAL CONDITIONS 

In this appendix the relations used in selecting initial conditions for the atom and 
molecule are derived. Present calculations involve collisions between molecules with 
specified internal energies and atoms selected at random from a heat bath. Conse- 
quently, initial conditions for the molecule involve only a random orientation of position 
and velocity vectors over a sphere, while initial conditions for the atom in the present 
collision model must be selected from a probability distribution corresponding to a 
Maxwell- Boltzmann flux distribution for the atoms at the specified heat bath tempera- 
ture. (The general method for choosing initial conditions from a known distribution is 
given in ref. 21. ) Implicit in the latter is the assumption that the translational tem- 
peratures of the molecules and atoms are the same. This situation corresponds to the 
common assumption in shock tube work that translational temperatures of all species 
are equalized in the shocked gas, while the internal degrees of freedom for molecules 
are in the process of "relaxing" toward equilibrium. 


Initial Conditions for the Atom 

An atom selected from the Maxwell- Boltzmann flux distribution with temperature T 
is placed on the surface of a sphere centered at the center of mass of the molecule, and 
with a radius equal to the interaction radius for the atom-molecule interaction (10 A). 
Without loss of generality the initial position was set at X^ = Xg = 0, Xg = -10 A (see 
fig. 1). Initial atom velocities were then chosen at random from 

N exp(-H A /kT)d£ (Cl) 

where N is a normalization factor to be determined. The flux distribution of atoms 
across the X^ - Xg plane is 


fdX = NL 


XeXP (" 2 ^ Ma > bcI 


1^1 


|xndx 


(C2) 


where B = I/kT and is is a unit vector along the X„ axis. One must determine N 

/ -*» /- -\ / r 3 
such that / fdX = 1 =sig •• X/ where \ig • X/is the average over the Maxwell- 

Boltzman flux distribution for X: 
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• ?>)^p(-^m AiBC |x| 2 )ck 


We introduce cylindrical coordinates for X by 


if j = |x|cos 0 

1*1 

Xg = X sin 0 > 
*3 = r 3 ' * 


(C3) 


(C4) 


The distribution f is then 

X 3 e X p[-i^ A BC (|x| 2 + JC 2 )]|x|d|x|d0dX 3 = jexp^Is 4 AjBC |x| 2 )|x|d|x|] 

X [exp^-i dK A|Bc x|x 3 dX 3 J(d0) (C5) 
where |x| goes from 0 to Letting |x| = s and a = /3ju^ gives 


1 

2 



exp(-as)ds = l/2a 


(C6) 


so that the normalized integral is 


2a 




= 1 


Also, 



-as)ds = -(l/2a)exp(-as 


)]q = (l/2a)(l - 



(C7) 
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Thus the normalized X distribution is 


\ ^A,BC^| 2 )l^l d l^l]|< 1 /^A,BC> ex pf \ ^BC^S ^3] 


x [d0/2ir] (C8) 


while the cumulative distributions are 


-\ 


(C9) 


Z i,2 =1 - ex p(-^^A,BC^ 12 


Z 3 ' 1 “ ex P^ ~ PV-A'BC *3) 


Z' Q = 0/2ir 


Setting the cumulative distribution functions (C9) equal to random numbers in the interval |_0, l] 
and recognizing that one minus a random number is also a random number in [0, l] gives 


0 = 2ttZ 


0 


i\ = (- 2 ln Z 1>2 /^ AjBC ) 1/2 


X3 - (-2 In z 3 /PM AjBC ) 


1/2 


(CIO) 


where Zq, Z^ ^ and Zg are random numbers. Thus, the initial conditions from 
(C4) become 


X, = (-2 in z 1: 2 /Pw A> Bc ) 1/2 cos( 2 itZ q ) 
X 2 = (-2 in Z h 2 /0«a, BC ) 1/2 S ln(2»Z 0 ) 

Xj = (-2 In ^/UMa.Bc)'^ 2 


(Cll) 
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Initial conditions are generated by selecting random numbers Z^ ^ Zg, and Zq from 
a linear distribution in the interval [0, 1] and substituting into equation (Cll). If this is 
repeated many times, initial kinetic energies H^(i) will be distributed according to a 
Maxwell-Boltzmann flux distribution at temperature T, as shown in table I. 


Initial Conditions for the Molecule, Given Hg^ and | Mg^, 


The orientation of Rg^ must be randomly distributed over a sphere. Thus, in- 
troducing the spherical polar coordinates 




Rg C i = |Rgc|sin e cos cp 


Rg C 2 = I R BC I sin e sin <P y 


R BC, 3 “ I r bcI cos 9 j 


and choosing |Rg C |(i)= |Rg C (e)|, we have 

cp = 27tZ 


<P 


9 = cos"*(l - 2Z 0 ) 


(C12) 


(C13) 


-q are random numbers in the interval [0,1], 


where Z „ and Z 

cp 

One can choose | RgQ I (i) anywhere between |Rg^ 


and R 


BC 1 


provided the 


max mm 

| r bc I distribution function can be integrated in closed form. If this cannot be done, 
approximations to the distribution must be made. Originally our program specified 
and H r and chose | Rg^, | (i) from a linear distribution between |Rjjq| and 

i — i , — , 2 min 

|R-d/-,I • However tabulation of results for | Mg^ | (i) = 0 trajectoris indicated no 


BC 1 


max 


detectable difference from tabulations where |Rg^| (i) = I R BC^ e ^ I * * n § enera l> the 
molecule underwent several complete vibrations before the atom came close enough to 
cause significant energy transfer, and it appears that statistical results are not sensi- 


tive to the choice of R 


BC 


|(i). 


For velocity components of the molecule, equation (Cl 2) is differentiated with respect 
to time, and with Vj^Rg^ l)(i) - 0, |5 b C | (i) is calculated from i m BC I R BC I 2 " H v 
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I 


I 


I III I 


I I I II I 1 1 ■■ III I 


I I 


III I 


where H y = H B( ~, - H r . A ± sign is then affixed to the result (by random number selec- 
tion) with equal probability of either choice. A plane of rotation is selected by choosing 

|e| = |4|z (ci4) 

1/2 

where Z is a random number and | 4 | = (2 H r / p BC | Rg C (e) | 2 ) is the magnitude of 
the instaneous angular velocity in the plane of rotation, and a ± sign is randomly affixed 
to |#|. Finally \cp\ is calculated from 

I M BC I /^BC = | R BC^I = 1^30(6)1 4 (^ 2 sin 2 0 + 0 2 ) 


or 


4 2 = (4 2 _ 0^)/sin 2 0 (C15) 

and a ± sign is randomly affixed to the result. With the initial components for 

-I- i -► | 0 oL/ 

and Rg£ thus chosen, the computer then calculates Hg C and |Mg C | and compares 
the results with the actual input before proceeding with the integration. In the unlikely 
event that | 0 1 =0 or j 4 | , and/or 0 = 0 or n, the computer rechooses 9, (p, and 9. 
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TABLE I. - FRACTIONAL POPULATIONS OF INITIAL ATOM ENERGIES FOR T = 1800 K 



. 01 
.02 
.03 
.04 
.05i 


.9476 
.8126 
.6538 
. 5102 
.4043 


.9359 
.8065 
.6616 
. 5302 
.4184 


.06 

.07 

.08 

.09 

.10 


.3058 

.2341 

.1731 

.1310 

.1002 


.3234 

.2426 

.1759 

.1268 

.0918 



- 


(t. 

(RCOL +NCOL )W^( RCOL +NCOL ) 

Required 

population 

W A 

/ j 

r aw 



Ju 


1 

1 


.9398 

.9378 


.8086 

.8071 


.6589 

.6600 


. 5240 

. 5218 


.4137 

.4027 


.3175 

.3053 


. 2397 

.2282 


.1749 

.1687 


.1282 

.1236 


.0947 

.0898 
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TABLE II. - RKR(J = 0) VIBRATIONAL ENERGY LEVELS 


FOR 79 Br 81 Br GROUND STATE (*£+) 


Vibrational quanta n- 

H v* D e 

Vibrational quantum 

H v* D e 

number. 

number, 


q 


q 


0 

0.01007 

40 

0.69485 

1 

.03011 

41 

.70839 

2 

.05000 

42 

.72170 

3 

.06974 

43 

. 73479 

4 

.08935 

44 

. 74764 

5 

. 10881 

45 

. 76026 

6 

. 12813 

46 

.77263 

7 

. 14731 

47 

. 78476 

8 

. 16635 

48 

. 79662 

9 

. 18525 

49 

.80823 

10 

.20401 

50 

.81957 

11 

.22263 

51 

.83063 

12 

.24111 

52 

.84142 

13 

.25944 

53 

.85191 

14 

.27763 

54 

.86212 

15 

.29568 

55 

.87202 

16 

.31359 

56 

.88162 

17 

.33135 

57 

.89090 

18 

.34896 

58 

.89986 

19 

.36643 

59 

.90849 

20 

.38375 

60 

.91679 

21 

.40091 

61 

.92474 

22 

.41792 

62 

.93233 

23 

.43478 

63 

.93957 

24 

.45147 

64 

.94643 

25 

.46801 

65 

.95292 

26 

.48439 

66 

.95903 

27 

. 50060 

67 

.96474 

28 

. 51664 

68 

.97004 

29 

. 53252 

69 

.97494 

30 

. 54821 

70 

.97941 

31 

. 56374 

71 

.98345 

32 

.57908 

72 

.98705 

33 

. 59423 

73 

.99020 

34 

.60920 

74 

.99289 

35 

.62398 

75 

.99512 

36 

.63856 

76 

.99686 

37 

.65294 

77 

.99811 

38 

.66712 

78 

.99886 

39 

.68109 





TABLE m. - PROBABILITY DISTRIBUTIONS FOR AH BC AND a|m bc | 2 PER COLLISION AT T = 1800 K 

[Entry n means given interval is within noise level; all entries except H 0C (i) = 0. 5 D g , |M BC | 2 (i) = 0 are based on 1000 RCOL; exception contains 1756 RCOL; estimate of reliability 
for fractional populations is ±0. 01. ] 


|M BC | 2 (i) £( RC0L > 

£(RCOL + NCOL) 


Fraction of (RCOL ry 

+ NCOL) dissociated W 1 = J f( AH BC )d(AH BC J 


w l/ w 2 



Integral lower limit, x 


ZJ 

- -30 

! 

_ 30 

- CO 


-n 

0. 0001 

0.001 

0.01 

0.1 

1.0 






Integral upper limit, 

y 






-1.0 

-0. 1 

-0.01 

-0.001 

-0.0001 

n 

GO 

00 

GO 

OO 

OO 


0 

0 

a 0. 308 

0 

A 

o/o 

0 /o 

0 /o 

0 /o 

0 /o 

0.76o/o;721 

0.24o/n 

0.189/0.279 

0.062/o.233 

0 /o. 173 

o/o.042 

0.1 

0 

0.318 

0 

l/l o/o 

0 /o 

0.000/0 

O.OOs/o 

0.015/0 

0. 758/). 713 

0. 227/n 

0. 176 /o. 287 

0.062/0.234 

0.00l/0.171 

0/0. 046 

0.1 

0.01 

0.315 

0 

A 

o/o 

0 /o 

0 /o 

O.OO8/0. 040 

0.028/n 

0.742/0.696 

0. 23o/n 

0.177/0.264 

0.06o/o.231 

0 /o. 169 

o/o.044 

0.1 

0.1 

0.315 

0 

l/l o/o 

0 /o 

0.001/0. 060 

0.024/0.076 

0.078/n 

0.712/0.688 

0.210/n 

0.167/0.235 

0.062/0.217 

O.OOo/o. 162 

o/o.039 

0.1 

1.0 

0.315 

0 

l 

o/o 

0 / 0. 096 

0.026/0. 137 

0.108/0. 144 

0. 139/n 

0.697/0.686 

0.164/n 

0.139/0.170 

0. 06l/o. 162 

0.000/0.130 

0/0. 048 

0.1 

2.0 

0.310 

0 

l/l o/o. 046 

0 /o. 125 

0. 066 /o. 156 

0.134/o. 162 

0.158/n 

0.699/0.691 

0. 143/n 

0.12l/0.147 

0.052/0.141 

0 /o. Ill 

o/o.037 

0.5 

0 

0.336 

0 

A 

7° 

0 /° 

O.Oll/o 

0.035/0 

0.052/0 

0.744/0.695 

0.204/n 

0.160/0.305 

0.058/0.240 

0 

CO 

0/0.058 

0.5 

0.01 

0.344 

0 

A 

o/o 

0 /o 

0.015/0 

0.038/0.042 

0, 065/n 

0.720/0. 669 

0.214/n 

0.162/0.289 

0.058/0.252 

0 /0.I86 

0/0. 056 

0.5 

0.1 

0.343 

0 

l/l o/o 

0 /o 

0.014/0. 066 

0.050/0. 088 

O.lOl/n 

0.692/0. 661 

0.207/n 

0.160/0.251 

0.06l/0.233 

0 /o. 173 

0/0.055 

0.5 

1.0 

0. 340 

0 

A 

7° 

0 /o. 099 

0.032/0.138 

O.lll/O. 145 

0. 15l/n 

0.669/0. 660 

0. 18o/n 

0.145/0.195 

0.068/0.189 

0 /0. 145 

0/0.059 

0.5 

2.0 

0. 337 

0 

V 1 

o/o. 048 

0 /o. 135 

0. O65/0. 171 

0.149/0.175 

0 .181 /n 

0. 671 /o. 664 

0. 148/n 

0.123/0. 161 

0.059/0.156 

0.000/0.123 

0/0. 050 

0.5 

5.0 

0.328 

0 

a 

o/o. 109 

O.OOl/o. 169 

0.12l/o.l98 

0. 177/0.202 

0.199/n 

0. 68o/o.673 

0.12l/n 

0.099/0.124 

0.04l/o.l20 

0 /0.095 

o/o.034 

0.5 

10.0 

0.315 

0 

l/l o/o. 118 

O.OO6/0. 175 

0. 131/0.200 

0. 180/0.205 

0.20l/i 

0.694/o.686 

0. 105/n 

0. 087/0.109 

0.042/o.l03 

0 /o. 085 

o/o.032 

0.9 

0 

0.409 

0 

V 1 °/° 

0 /° 

0.036/o 

0.093/0 

0.138/o 

0.655/0.617 

0.208/n 

0, 160/0. 383 

0.072/0.293 

0 

CO 

0/0. 089 

0.9 

0.01 

0. 399 

0 

V 1 °/° 

0 /° 

0.038/0 

0.095/0.055 

0, 15o/n 

0,654/0.612 

0.196/n 

0.145/0.332 

0.06l/0.286 

0 /o. 206 

0/0. 088 

0.9 

0.1 

0.397 

0 

l/l o/o 

O.OOO/O 

0.039/0,080 

0.099/0.098 

0.169/n 

0.625/0. 609 

0.205/n 

0.144/0.292 

0.062/0.277 

0 /0.203 

0/0. 085 

0.9 

1.0 

0.392 

0 

A 

o/o 

O.OOl/o. Ill 

0.044/0.158 

0. m/o. 165 

0.188/n 

0.617/0.610 

0. 195/n 

0.150/0.225 

0.07l/0.221 

0 /o. 173 

0 /o. 087 

0.9 

2.0 

0.384 

0 

l/l o/o. 051 

0,00l/0.131 

O.O66/0. 172 

0. 145/0. 177 

0.187/n 

0.627/0. 616 

0.187/n 

0.152/0. 206 

0,069/0.202 

0 /o. 163 

0/0. 081 

0.9 

5.0 

0.368 

0 

A 

o/o, no 

0,002/0. 180 

0.117/0.207 

0.187/0.213 

0.215/n 

0. 643/o.632 

0.142/n 

0.119/0.154 

0.055/o.l49 

0.000/).125 

0/0. 063 

0.9 

' 10.0 

0.340 

0 

A 

o/o. 148 

O.Ols/o, 196 

0.153/0. 216 

0.204/0. 221 

0.227/n 

0. 67o/o,661 

0. 103/n 

0.08l/0.118 

0.036/0.113 

0 /0.O88 

0/0.039 

0.999 

0 

a 0. 658 

0.316 

l/l o/o 

0 /o 

0 /o 

O.OOl/o 

0.055/0 

| 0.498/o. 514 

0.448/n 

0.122/0.486 

0.028/o.324 

0 /0. 186 

0/0. 082 

1.000 

' 0.1 
1 

a 0. 658 

0.315 

A 

o/o 

0 /o 

0 /o. 151 

0.003/0.220 

0. 08o/n 

0.482/0. 362 

0.438/n 

0.138/0.418 

0.034/0.352 

0 /0.211 

0/0.095 

1.009 

1.0 

a 0. 658 

0.515 

l/l o/o 

0 /0. 120 

0 /o.266 

0.034/0.292 

j 0.178/n 

0.493/0.350 

0.329/n 

0.128/0.357 

0. 04o/o. 327 

0 /o. 209 

0/0. 086 

1.019 

2.0 

a 0, 658 

0.608 

A 

o/o. 016 

0 /o. 140 

0 /o. 287 

0.044/0.322 

0.210/n 

0.53l/0.351 

0.259/n 

0.10l/0.327 

0.038/0.290 

0 /o. 169 

0/0. 081 

1.049 

5.0 

a 0. 658 

0.650 

iA 

o/o. 010 

0 /0.057 

0.002/0.257 

0.015/0.340 

0.145/n; 

0.676/0.368 

0.178/n i 

0.038/o. 292 

0.024/0.212 

0 / 0. 063 

o/o.032 

1.099 | 

10.0 

a 0. 658 

0.633 

A 

0/0.034 

0 /o. 133 

0.020/0. 289 j 

0.060/0.33^ 

0.246/n 

0.512/0.356 

0. 242/n 

0.075/0.313 

0.04l/o. 277 

O.OOl/o.l53 

o/o.053 


CO 

-3 


Trajectories determined by fixed step size. 



TABLE IV, - PROBABILITY DISTRIBUTIONS FOR AH BC AND a|M bc | 2 PER COLLISION THROUGH A RANGE OF HEAT BATH TEMPERATURES 

1 — i o 

Mpcl (i) 0. T(A r ) - 1800 K are based on 1000 RCOL; an estimate of reliability 

for fractional populations is ±0.01.] 


H BC (i) ||M R J 2 (i) 


£(RCOL) 


' DRCOL + NCOL) 


Temperature, 

K 


j: 


■ fj <-i«bcI 2 >(-I“bcI 2 ) 


Integral lower limit, y 







-- 

-- 


-- 

-n 

0.0001 

0.001 

0.01 

0.1 

1.0 

10.0 








Integral upper limit, y 






-1.0 

-0.1 

-0.01 

-0.001 ■ 

-0.0001 

n 


- 

! 

QO 

00 


oO 

0 

0 

a 0. 309 

300 

o/o 

0 /o 

0 jo 

0 /o 

0 A 

0.719/o.697 

! 0. 28l/n 

1 0. no jo. 303 

0.007/0.272 

0 jo. 140 

o/o.OOl 

0/0 

0 

0 

a 0. 309 

1000 

ow— 

0 

0 

1 0 jo 

0 / 

0 


^ jo 

0. 744/0.709 

0.256/n 

0. 193/0.291 

0. 042/0. 249 

"5 /0 . 173 

0/0.019 

ojo 

0 

0 

a 0. 308 

1800 

0/0 

0 

/° 

0 / 0 

0 / 

0 


0 jo 

0. 760/o. 721 

0. 24o/n 

0. 189/0.279 

0.062/0.233 

0 jo. 173 

0/0 042 

0/0 

0 

0 

a 0. 309 

3000 

o/o 

0 


0 jo 

0 / 

r~ 

0 


0 jo 

0.77o/o.732 

0.23o/n 

0. 186/o.268 

0.083/0.223 

O.OO2/0. 172 

0/0. 058 

ojo 

0 

0 

a 0. 309 

5000 

o/o 

0 


0 /o 

0 / 

0 


0 /o 

0. 779/0.747 

0.22l/n 

0. 184/0.253 

0.099/0.215 

o.ooe/o.172 

0/0. 076 

0/0.001 

0 

0 

a 0. 309 

10000 

o/o 

0 

b 

0 jo 

0 / 

1 

0 


0 jo 

0.787/o. 762 

0.213/n 

0.179/0.238 

0.117/0.205 

0.024/o. 168 

o/o.097 

0/0. 008 

0. 5 

0 

0.340 

300 

o/o 

0 

b 

0 b 

0.018/ 

T 


0.044/0 

0.697/0.668 

0.258/n 

0. 162/0. 332 

0.007/0.296 

0 /o. 171 

0/0. 006 

ojo 

0.5 

0 

0.347 

1000 

o/o 

0 

b 

0.004/0 

0.028/ 

0 


0.047/0 

0. 716/0.671 

0. 236/n 

0.177/0.329 

0.040/0. 272 

0 /o. 190 

0/0.035 

ojo 

0. 5 

0 

0.336 

1800 

o/o 

0 

b 

0.011/0 

0.035/ 

^0 


0.052/0 

0.744/0.695 

0. 204/n 

0. 160/0. 305 

0.058/o. 240 

0 / 0.181 

o/o.058 

ojo 

0.5 

0 

0.349 

3000 

o/o 

0 

b 

0.018/0 

0. 048 j 

io 


0.072/0 

0.737/0.695 

0. 19l/n 

0. 156 /o. 305 

0.075/0.237 

O.OO3/0. 186 

o/o.077 

0/0.001 

0.5 

0 

0.349 

10000 

o/o 

0. 

000/0 

0.026/0 

0.056/ 

^0 


0.102/0 

0.72l/o.735 

0. 177/n 

0.152/0.265 

0.108/0.217 

0.023/0.179 

o/o.l07 

0/0.014 

0.5 

1.0 

0.345 

300 

o/o 

0 

/o. 188 

0.044/0.213 

0.198y 

lo. 

216 

0.219/n 

0. 660/0. 656 

0. 12l/n 

0.093/0.128 

0.004/0.123 

0 jo. 092 

o/o.005 

0/0 

0.5 

1.0 

0.344 

1000 

o/o 

0 

Jo.xzi 

0.037/0.170 

0. 14 ly 

/ST 

174 

0.174/n 

0.664/0.657 

0. 163/n 

0. 124/0.170 

0.044/0. 165 

0 /o. 129 

0/0. 038 

ojo 

0.5 

1.0 

0.340 

1800 

o/o 

0 

1 0.099 

0.032/0.138 

0 . Illy 

lo. 145 

0. 15l/n 

O.669/0. 660 

0. 18o/n 

0. 145/o. 195 

0.068/o.l89 

0 /0.145 

o/o.059 

ojo 

0.5 

1.0 

0.346 

3000 

o/o 

0 

/ 0.089 

0.032/0.143 

0. 103^ 

/o. 151 

0.149/n 

0. 671 jo. 657 

0. 1 80/ n 

0. 147/o. 192 

0.079/0. 187 

0.004/0.147 

0/0. 076 

0/0.001 

0.5 

1.0 

0.346 

10000 

o/o 

0.002/0.056 

0.032/0.105 

0.074^ 


121 

0.129/n 

0.679/o. 658 

0.192/n 

0.155/0.221 

O.lll/o. 206 

0.026/0. 155 

o/o.l03 

0/0. 018 

0.5 

5.0 

0.328 

1800 

o/o. 109 

0. 

001/0. 169 

0.121/0. 198 

0.177 ( 

K 

202 

0.199/n 

0.680/0.673 

0 .12l/n 

0. 099/0.124 

0.04l/0.120 

0 /o. 095 

o/o.034 

0/0 

0.5 

5.0 

0.328 

10000 

o/o. 060 

0. 

003/0. 102 

0.O63/0.142 

0. 104y 

lo. 

149 

0.139/n 

0.685/0.673 

0 . 176/n 

0. 142/o. 178 

O.lOl/o. 171 

0.028/o. 133 

0/0. 089 

0/0.019 

0. 5 

10.0 

0.315 

1800 

o/o. 118 

0. 

006/0. 175 

0.131/0.200 

0.180^ 

Jo. 

205 

0.20l/n 

0.694/0.686 

0. 105/n 

O.O87/0. 109 

0.042/o. 103 

0 jo. 085 

o/o.032 

ojo 

0.5 

10.0 

0.314 

10000 

0 jo. 086 

0. 

009/0. 123 

0.091/0.157 

0.125 j 

fa- 

163 

0. 154/n 

0.695/0.686 

0. 15l/n 

0.125/0.151 

0.092/o.l45 

0.024/o.U5 

o/o.079 

0/0.013 

0.5 

0.25 

0.343 

300 

o/o 

0 

jd. 078 

0.001/0.139 

0.109^ 

fa. 

146 

0. 156 /n 

O.668/0. 658 

0. 176/n 

0.128/0. 196 

0,007/0.188 

0 / 0.133 

0/0.004 

ojo 

0.5 

0.82 

0. 343 

1000 

o/o 

0 

/o. 115 

0.029/0.152 

0.135 

/o. 158 

0.166/n 

0.669/0.659 

0. 16s/n 

0.132/o. 183 

0.042/o.l75 

0 / 0.135 

o/o.037 

0/0 

0. 5 

1.48 

0,341 

1800 

o/o. 031 

0 

jo. 123 

0.049/0.159 

0.134 

jo. 165 

1 0.166/n 

0.669/0. 660 

0. 165/n 

0. 134/o. 175 

1 0.063/0. 169 

O.OOo/o.l33 

o/o.057 

0/0 

0.5 

8.24 

0.321 

10000 

o/o. 082 

! 0. 

OO7/0. 119 

0.082/0.152 

0.121 
t 

E 

161 

0. 15l/n 

0. 69o/o. 680 

0. 159/n 

0.13l/o.l59 0.092/0.152 0.025/0.120 o/o.082 

1 0/0.015 


Trajectories determined by fixed step size. 



TABLE V. 


CORRESPONDENCE (FRACTION OF RCOL) BETWEEN GAINS AND 


LOSSES IN H 


BC 


AND M 




FOR T = 1800 K 


^BCI 5 


CQ 

5 



Gains 

Losses 

CQ 

C 



rt 

O 

(g,g) 

(g,D 

CQ 



CQ 

CQ 

o 

a,g) 

d,i) 


g = Gains 
l = Losses 


( h bc^’ I^bcI 2 ^) 


1.000 

0 

0 

0 

(0,0) 

0. 849 

0 

0. 151 

0 

(0.1 D e ,0) 

0.762 

0 

0.238 

0 

(0.5 D e ,0) 

0.353 

0.030 

0.029 

0.588 

(0.5 D e , 5.0) 

0.433 

0.075 

0. 144 

0. 348 

(0.9 D e , 1.0) 

0.548 

0. 190 

0.104 

0. 158 


(1.000 D e ,0. 1) 


0. 782 

0.011 

0. 075 

0. 132 

(0.1 D e ,0.01) 

0.685 

0.023 

0.172 

0.120 

(0.5 D e , 0.01) 

0.320 

0. 022 

0.027 

0.631 

(0.5 D e , 10.0) 

0.412 

0.088 

0. 126 

0. 374 

(0.9 D e , 2.0) 

0.473 

0.121 

0.078 

0.328 


(1.009 D e , 1.0) 


0. 700 

0. Oil 

0. 051 

0.238 

(0.1 D e ,0.1) 

0,612 

0.039 

0.125 

0.224 

(0.5D e ,0.1) 

0. 593 

0 

0.407 

0 

(0.9D e ,0) 

0.306 

0.094 

0. 114 

0.486 

(0.9 D e ,5.0) 

0.419 

0.080 

0.089 

0.412 


(1.019 D e , 2.0) 


0. 528 

0. 017 

0. 015 

0. 440 

(0.1 D e , 1.0) 

0, 507 

0.036 

0,068 

0. 389 

(0.5D e , 1.0) 

0. 504 

0.045 

0.344 

0.107 

(0.9 D e , 0.01) 

0. 261 

0.054 

0.087 

0. 598 

(0.9D e , 10.0) 

0. 338 

0.093 

0.130 

0.439 


(1.049 D e 5.0) 


0.470 

0.004 

0.007 

0. 519 

(0.1D e ,2.0) 

0.420 

0.035 

0.059 

0. 486 

(0.5D e ,2.0) 

0.485 

0.055 

0.262 

0. 198 

(0.9D e ,0.1) 

0.790 

0 

0.210 

0 

(0.999 D e ,0) 

0.428 

0.037 

0.056 

0.479 


(1.099 D e , 10.0) 
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TABLE VI. - CORRESPONDENCE (FRACTION OF RCOL) BETWEEN LARGE AND SMALL GAINS AND LOSSES 


IN H 0C AND |m bc | 2 FOR T= 1800 K 



0 I 0 0 1 0 


( 0 , 0 ) 




(O.5D 0 ,O) (0.5D e ,0.01) (0.5D 0.1) (0. 5 D 1.0) (0.5D 2.0) 














0.085 jo. 150 

0.014 jo. 010 

1 

0 [0.118 

| 

i 

0.001 jo. 005 

1 

0.006 10.014 
1 

0.086 jO. 141 

i 

0.001 |0. 008 

I 

0.002 |0. 359 

(0.5D ( 

5.0) 

( 

0.081 lo. 186 
__L 

1 

0. 021 J 0. 040 

1 

0.007 [0. 159 

1 

0.003 [0.011 
1 

0.041 lo. 051 

1 

0. 112 0.176 

i 

1 

0.007 jO. 045, 
i ! 

1 

0.004 |0. 056 
1 

(0.9 D 

e .i.°) 

“i 

0.240 Jo. 257 

1 

0. 190 1 0 

j 

0 0.051 

0 1 0 
1 

0.092 jo. 012 

0. 158 I 0 
| 

1 

0 1 0 
1 

0 j 0 


0.063 jo. 128 
1 

0.001 jo. 128 

0.010j0.009 

1 

0.00110.002 

1 

1 

0.01410.013 

1 

0.086] 0.130 
1 

1 

0 1 0 
] 

1 

0 1 0. 415 

1 

(0. 5 D £ 

.,10.0) 

1 1 
0. 070 jo. 182 0.022 jo. 047 

i i 

0.002 |0. 158 0.004(0.015 | 
1 1 

1 

0.038 |0.045 
1 

1 

0. 092 | 0. 153 
1 

1 

0.005)0.038 
1 

0. 004 jo. 125 
1 


(0.9 D e , 2.0} 


1 

0. 16510.247 0. 107 1 
i 

0.014 

CD 

O 

o 

° 

0 

0 

0. 069)0.009 

0.160 

jo. 168 

0 j 0 
1 

0 

! ° 

i 


(1.000 D .0.1) 


(1.009 D e> 1.0) 


i 

0. 238 jo. 179 

0 


1 

0.007 JO. 169 
1 

0 

0 



0.220 [0.098 
1 



1 

0.011 10.078 
1 




(0.9 D e ,0) 


1 1 
0.049 |0.128 0.026 '0.047 


1 

03 

eg 

O 

o 

0.001 10.020 
l 

l 

0.026 [0.050 
1 

0.060 Jo. 146 

J 

0.004 0.034 
1 

0.004)0.276 
1 


(0.9 D ,5.0) 


0 . 169 lo. 192 0. 06310.017 
I I 


! 

0 [0.058 

j 

1 

0 I 0 

0.083 jo. 006 

0.216 0.196 

1 

1 

0 j 0 


0 | 0 
1 


~~ 1 

0. 173*0. 179 

r 

0. 043 j 

0 

1 

0. 005 j 0. 147 

1 

0.002) 

1 

0 

0. 142| 0. 108 
1 

0.105| 
l_ 

1 0 

J 

0.011|0.083 
1 

© 

0 

0 

CO 

0 


(0.9 D e ,0.01) 


1 

0.064 jo. 103 

1 

0.015I0.028I 
1 1 

1 

0.002 |0.092 
1 

0.001 

0.010 

0. 021 |o. 047 
1 

0.058 

0.109 

0. 001 jo. 018 

1 

0.001 

0.430 

(0.9D e 

,10.0) 


1 

0.24710.054 

1 

0.088 

0.005 

1 

0 10.037 

0 

0 

0. 126 0.004 
1 

0.357 

0.079 

0 j 0 
1 

0 

jo. 003 
| 


1 

0. 143 [ 0 . 190 

r 

0.050 j 

0 

0.005 lo. 147 

l 

1 

0.005 | 
1 

0 

1 

0.080 |0.092 
1 

1 

0. 191 1 

0 

1 

0.008 jO. 082 
1 

f 

0.007 | 

L 

0 


(0.9D e ,0.1) 


1 

0.515 10.232 
1 

1 

0 1 0 
1 

1 

0 jO. 043 

1 

0 I 0 
1 

1 

0.202 |0. 008 
1 

0 0 

0 j 0 

1 

0 j 0 

_j 

(0.999 D e> 0) 

1 

0. 196 lO. 169 

I 

0.035 |0.002 
1 

1 

0 |0.063 

1 

0 j 0 

1 

1 

0.055)0.001 

1 

0.279 10.170 

j 

I 

— 

0 1 0 
1 

0 ( 0.030 

1 


(1.019 D e> 2.0) 


(1.049 D ,5.0) 


(1.099 D ,10.0) 











TABLE VH. - CORRESPONDENCE (FRACTION OF RCOL) BETWEEN LARGE AND SMALL GAINS AND LOSSES 



(0,0,10 000) 


(0. 5 D g) 0,300) 


(0. 5 D e , 0,1000) 


(0. 5 D e , 0,1800) 


(0. 5 D e> 0,3000) 























































TABLE VIII. - VALUES OF |M BC | 2 (i) REQUIRED FOR AHg C AND a|m bc | 2 GAIN-LOSS 


BALANCES AND COMPARISON WITH 


<Kc! 2 > 


FOR EQUILIBRIUM RIGID ROTATORS 


H BC (l) 

Temperature, 

^BC 

(i) for 

<^|MBc|^f° r rigid rotators with 

I r bcP ’ 


T, 

vibrating rotator 

I^BC^I 

a 5. 27(1) 

a 5. 62(5) 

a 8. 9(7) 

/Ir l 2 \ 


K 

Balance for, 

ah bc 

Balance for, 

a K c I 2 

\I R BcI /min 

0.1 D e 

1 800 

1. 53±0. 06 

1. 50±0. 06 

0.286 

0.290 

0.309 

0.493 

0.238 

• 5 D e 

300 

. 30±0. 03 

. 36±0. 06 

.477 

.483 

. 515 

.821 

.397 

• 5D e 

1 000 

. 79i0. 02 

. 98±0. 02 

.954 

.966 

o 

CO 

1.64 

. 794 

• 5D e 

1 800 

1. 45±0. 02 

1. 73±0. 02 

1.72 

1. 74 

1.86 

2.96 

1.43 

.5D e 

10 000 

9. 46±0. 03 

7. 90±0. 03 

4.77 

4.83 

5. 15 

8.21 

3.97 

• 9 D e 

1 800 

2. 00±0. 00 

2. 30±0. 11 

9.54 

9.66 

10.3 

16.4 

7.94 


Rotational 

temperature, 

K 


300 
500 
1 000 
1 800 
5 000 
10 000 


a For the vibrating- rotator /IR q^I 2 ^ = 5. 27(1) A 2 at 0. 1 D_, 5.62(5) P? at 0. 5 D , 8.9(7)A‘ 

cussion accompanying eq. (9)). ( \ ~ 4.33 A2 calculated from average of RgM at 

\| / mir, ] Imin 


at 0. 9 D e (see dis- 


and 0. 5 D e . 
V 1 = 0.5D e: 


Because of vibration- rotation coupling, 
' 2 0. 


except for M 1 


BC 


at H 


BC 


0.5 D^ 


mm 

is not equal to 


V x =0.1D e 


BC 


at 


BC | 
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TABLE IX. - PROBABILITY DISTRIBUTIONS AND MATRIX TABULATIONS OF AH BC AND a|m bc | 2 FOR TRAJECTORIES CALCULATED WITH ATTRACTIVE PART OF Vg REMOVED 

* 2- fj 



Integral lower limit, x 


- or> 

i 

- | - 


-n 

0.0001 

0.001 

0.01 

0.1 

1.0 

10.0 









Integral upper limit, y 










-1.0 

-0.1 -0.01 

-0.001 

-0.0001 

n 

- 

CO 

X) 

- 

- 

QO 

0.5 

0.25 

0.348 

300 

o/o 

0 /0.062 0 /. 150 

0. 07/6. 170 

0. 17l/n 

0.69/6,659 

0.13/n 

0. 076/0. 171 

0.00/6.152 

0 /. 070 

0/6. 001 

0/0 

0. 5 

1.48 

0.342 

1800 

: 0/0.021 

0 /.\20 0.03/0.166 

0. 132/0. 174 

0. 171/n 

0. 67S/o. 660 

0. 154/n 

0. l 20 /o. 166 

0.04/0.155 

0.000/6. 116 

O/O. 040 

0/0 

0.5 

8.24 

0.319 

10000 

o/6. 075 

0. 00 f/. 129 j 0. 079/0. 154 

0.130/0.158 

0. 15l/n 

0.690/0.682 

o.is/n 

0.141/0.159 

0.08/0.154 

0.021/6.131 

/0.080 

o/o.012 


BCl 


(0 


£(RCOL) Temperature, 

E(RCOL+NCOL) T ’ 

K 


vf 


1 

0. 251 10. 195 

°_j 0 

0 0.007 

1 

0 j 0 

1 

0. 051 1 0 

1 

0. 318 1 0. 718 


, 1_ , 

1 -j 1 


| 0 j 0 

, 0 1 0 
1 1 1 


(0. 5D e> 0.25, 300) 


0. 132 1 0 . 188 

0.009 0.011 

i 


0 jo. 137 

0 10.001 

1 

1 

0. 017 1 0. 010 
1 

1 

0.150(0.247 

1 

10.001 10.004 

- 

1 

0 |0. 093 1 


(0. 5 D ,1.48, 1800) 


0.08010.129 0.013 0.014 


0.003 jo . 266 

_ 1 

0 

0.010 

1 

0.008 10.010 
t | 

0.081 

0 . 139 

| 0 lo . 004 

0.001 

0 . 242 


(0. 5 D e ,8.24, 10 000) 


cn 












Dissociation events with ljtf BC l 2 loss 
Dissociation events with |^fgc| 2 9 a * 


Fraction of 
total events 



-n 0 n 
AH bc or A|M B q| 2 


Figure 6. - Schematic drawing of a family of probability distributions for 
internal energy changes, AH BC , or angular momentum changes, 
A|M B c| 2 - P er collision. 



|/V?Bc| 2(i) 


Figure 7. - Comparison of dissociation events with gain or loss in (T * 1800 K). 
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Fraction of total events 


«bc<». |M B c| 2(i,) 


Losses 


Gains 


T(Ar) 


(0. 5 De, 1.0) 


(0.5 De.O) 


(0.5 De, 1.0) 



T(Ar) * (0 * 5De - 0) 


(0. 5 De, 1.0) 


T(Ar) 


J0JDe,0)| 


(H Bd^ l^Bc! 2<f)) 


"( 0 , 0 ) 


ZZ^-rS 6 ’?) Noise level 
(0. 5 De, 1.0) 



. T(Ar) 

( 0 , 0 ) 

(0. 5 De, 0) 0.0001 to «> 

■ T(Ar) 

-( 0 , 0 ) 


— 10.5De,0) 

(0.5De, 1.0) 0.001 to °° 


( 0 , 0 ) 

(0. 5 De, 0) 
(0. 5 De. ; ' 


■ T(Ar) 



0. 01 to 00 


0. 1 to°o 


T(Ar) ■ 


T(Ar) 


(a) Behavior of AH BC probability distributions as function of heat bath temperature. 
Figure 8. - Temperature dependence of probability distributions. 
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Figure 8. - Concluded. 
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